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We show that an ensemble of identical d-level atoms can be efficiently described by d — 1 col- 
lective oscillator degrees of freedom in the vicinity of a product state with all atoms in the same, 
but otherwise arbitrary single-particle state. We apply our description to two different kinds of 
spin squeezing: (i) when each spin-f atom is individually squeezed without creating interatomic 
entanglement and (ii) when a particular collective atomic oscillator mode is squeezed via quantum 
non-demolition (QND) measurement and feedback. When combined in sequence, the order of the 
two methods is relevant in the final degree of squeezing. We also discuss the role of the two kinds 
of squeezing when multi-sublevel atoms are used as quantum memories for light. 
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I. INTRODUCTION 

Large ensembles of identical particles are excellent can- 
didates for light-matter interfaces, quantum memories, 
repeaters for long distance quantum communication, and 
registers for quantum computing. This is not only due to 
the fact that in these systems certain collective quantum 
degrees of freedom efficiently couple to the electromag- 
netic radiation field, but also because the desired control 
is provided by simple, experimentally accessible interac- 
tion mechanisms. This is the case, e.g., for large collec- 
tions of N identical spin-^ particles. Such an ensemble 
can be effectively described using a single collective spin, 
if the initial state is a pure product state with all parti- 
cles prepared in the same single-particle state, and if the 
Hamiltonian describing the system is a sum of identical 
single-particle operators. Since, in this case, the quan- 
tum state of the entire collection remains invariant under 
permutations of particles, the corresponding restricted 
Hilbert subspace is equivalent to the states of a large 
spin-^ spin. If, in addition, the dynamics only weakly 
perturbs the initial state, the collective states explore 
only a limited range of states, and the pseudospin eigen- 
basis can be mapped onto an oscillator basis of states, 
and one may benefit from the simple and well studied 
properties of harmonic oscillator systems. 

In this paper we wish to generalize the oscillator de- 
scription of two-level systems to multilevel particles. A 
natural example is the one of atoms having a Zeeman 
degenerate ground state with total angular momentum 

or a hyperfine ground state manifold with a range of 
angular momentum quantum numbers F' and associated 
\M'\ < F' . In a large ensemble of such atoms, the col- 
lective spin picture is less useful, and we will proceed 
directly to an effective oscillator description of the sym- 
metric collective states of the system and of the system 
dynamics. Various methods allow the control of single 
particle states, and the atoms in a large ensemble can be 
prepared in essentially any superposition of the hyperfine 
ground states Pl. Therefore, it is a particular purpose of 



our work to develop a theory of collective states in the 
vicinity of a product state of arbitrary single atom states, 
and to determine the interplay between the collective and 
single particle properties of such samples. 

In Sec. |TT] we review the usual Holstein-Primakoff ap- 
proximation, describing a collection of spin-^ particles 
by a single harmonic oscillator degree of freedom, and 
we generalize this description to multilevel systems ex- 
panded around product states of arbitrary single particle 
state vectors. We discuss the representation of collective 
operators, which are a sum of single particle operators 
over the entire ensemble, and we demonstrate that in the 
vicinity of product states, such operators can be associ- 
ated with oscillator quadrature operators in a generalized 
Holstein-Primakoff approximation. In Sec. Ill we turn 
to a special discussion of spin squeezing, and we point out 
the formal distinction between the effect of spin squeezing 
within each F > ^ hyperfine angular momentum mani- 
fold and squeezing of collective spin degrees of freedom. 
We also show that, when applied sequentially, the order 
of the two kinds of squeezing is relevant. In Sec. |IV| 
we discuss the use of samples of multi-sublevel atoms as 
quantum memories for light and the role of internal spin 
squeezing in such atoms. Sec. |V] concludes the paper. 



II. MULTILEVEL HOLSTEIN-PRIMAKOFF 
APPROXIMATION 



A convenient way to describe an ensemble of spin-^ 
particles or two-state atoms which are homogeneously 
coupled to external perturbations is in terms of the col- 
lective spin operator J = X^jLi s*--'-' • If the initial state 
of the system is a product state with all members oc- 
cupying the same pure state, the dynamically accessible 
Hilbert space is characterized by a single angular mo- 
mentum ladder. It is then straightforward to introduce 
oscillator-like creation and annihilation operators in the 
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Holstein-PrimakofF representation, 

Jx- = a} \/2J — a^a, = J — a) a, (1) 

where J is the total angular momentum quantum num- 
ber, which is a constant of motion in the cases considered, 
and the x quantization axis is chosen to be the direction 
of polarization in the initial product state. For macro- 
scopic polarization {J 1) and if the system stays in the 
vicinity of the spin coherent state | J, m j = J) , the atomic 
oscillator picture is especially efficient and it allows us to 
directly define the quasi-canonical atomic quadrature op- 
erators, 



X = 



1 



/J 



V' 



P = 



1 



/J 



(2) 



with [X, P] = iJx/J ~ i- 

For an ensemble consisting of atoms with larger angu- 
lar momenta F > ^ , it is also possible to define the collec- 
tive spin operator J — ^2^=1 ^^"'^ ■ 1^ ^^^h atom is initially 
prepared in the same spin coherent state, and the Hamil- 
ton operator is an element of the operator algebra gen- 
erated by the components of J, this collective spin again 
provides a general description of the system. In that case, 
it is possible to think of each atom as a collection of 2F 
fictitious spin-i particles, and to think of the initial state 
as a symmetric product state of the total of 2FN ficti- 
tious spin-1 particles, which evolves in the manifold of 
states of the large collective spin with J = NF. The col- 
lective spin components constituting the Hamiltonian are 
invariant under the exchange of any two fictitious spin- 1 
particles irrespective of whether they belong to the same 
atom or to different atoms. The collective spin descrip- 
tion is, however, incomplete in general. For example, 
the atoms may be individually prepared in arbitrary in- 
ternal superposition states J2m Cml-F, rnp) corresponding 
to correlations (entanglement) among their own fictitious 
spin-1 constituents. Even if the atoms are thus all pre- 
pared in the same state, the 2FN fictitious spins are not 
equivalent: those pertaining to the same atom share cor- 
relations in the initial atomic product state, while the 
ones that pertain to different atoms do not. 

We will now introduce a generalization of the Holstein- 
Primakoff representation that is capable of treating not 
only the total angular momentum operators but any per- 
mutation invariant sum of single-particle operators of an 
ensemble of d- level systems. 



A. Generalized spin operators and collective 
atomic oscillators 

Let us consider an ensemble of N identical d- level sys- 
tems, and let {|'/'Q)}a=o denote an arbitrary orthonormal 
single-particle basis. Throughout the paper, we restrict 
ourselves to symmetric collective states, i.e., those in- 
variant under permutation of the internal state of any 
two particles. For spin-^ particles, this corresponds to 



the maximal total angular momentum manifold, and our 
representation coincides with the standard approach. 

In the general case, the symmetric subspace is spanned 
by the occupation number states 



|no,mi,;2,---) = 



1 



V n\ m! U . 

X \4>l)n+l ■ ■ ■ \4>l)n+m \4'2)n+m+l ' ' 



perm 



|0o)l • • • \4'o)n 



(3) 



which means that n atoms are in the internal state \(j)o), 
m in \4>i), I in |(/>2), etc., and they sum up to iV = n -I- 
m + I + . . . atoms. It is easy to see that the collective 
operators 



N 

Zap = ^ |'^a)ji((/'/3|, 



(«,/? = 0,1,. 



1) (4) 



keep the symmetry of the sub-Hilbert space in consider- 
ation. Furthermore, we observe the following properties, 



T^apWa^mp, . . .) = V(n + l)m 

X |(n+ l)„,(m- 1)^,...), (5) 
ilaaWa.mp,...) ^n\na,mp,...), (6) 
and the commutator identity, 

[S(j-y, S^^] — Sq,^ — (5q.^S^^. (7) 

All other commutators, which can not be brought on the 
form [Sq,^,S^^], vanish. 

Any pure state of a spin-i particle is a spin coherent 
state, i.e., a spin up state along a suitably chosen axis, 
say X. The corresponding ensemble state |J, TOj = J) 
serves as the natural reference state: it is the ground 
state of the effective atomic oscillator in the Holstein- 
Primakoff approximation. With a general single-atom 
basis vector |0o), we take the product state 0^ |(/)o) as 
the reference state and introduce d — 1 independent col- 
lective atomic oscillator modes with creation and annihi- 
lation operators that redistribute the atomic populations 
between |(/)o) and the other basis states \4>a), 

al|no,mQ, . . .) = ^/m + l \{n - l)o, (m + 1)^, • • •), (8) 
aJjadncmQ, . . .) ^ to |no,mQ, . . .), (9) 

with a = 1, . . . ,d — 1. The reference state itself corre- 
sponds to the multi-mode vacuum state of the oscillators, 
while the excitation number states coincide with the sym- 
metric states ^ . In analogy with the Holstein-Primakoff 
representation, directly comparing Eq. ([S]) with (|8| and 
Eq. ^ with (|9|, respectively, the generalized collective 
spin operators can be expressed as follows. 



Sa/3 = aLa/3, ^ao = N -J2p^o^lo-i3, (10) 
furthermore, using Eq. 0, 



■^00 



N ■ 



J2dida. (11) 
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B. Linearization around the reference state: 
Holstein-Primakoff expansion of collective operators 

In the vicinity of the reference state |0o), i-e., for 
small number of excitations of the atomic oscillators, only 
the population Eqo has macroscopic expectation value, 
and the spin operators can be approximated as 

Sa/3 = a^ap < ±M-VNai < Eqo ~ N. (12) 

This allows us to simplify collective operators that are 
permutation invariant sums of single-particle operators, 
or belong to the algebra generated by such operators. 

Since {\4>a){(f'i3\}af3 forms an orthonormal single- 
particle operator basis, for any Hermitian single-particle 
operator O^^'' the corresponding iV-body operator can be 
expanded as 



N 



3 = 1 ap 



(13) 



where we used the shorthand Oai3 ^ {(j)a\0^^^\(j)p). Now 
we distinguish between two cases. If |(/)o) is an eigenstate 
of 0'^^\ then OaQ = for all a ^ and the collective 
operator becomes a c-number plus a correction that is 
quadratic in the creation and annihilation operators. 



00 



^ {O^a - O^^S^pYaiap. (14) 

a/3^0 



For example, taking the spin projection eigenstates 
\4>a) = \F,m^p ~ F — a), the longitudinal component of 
the total angular momentum is Jx = X]a(-^ ~ o:)'Eaa = 
Ea^^oa^aa- 
If |0o) is not an eigenstate of O*^^'', we get cross terms 
between the a — and the a ^ components, and we 
obtain the dominant contribution 



O « NO. 



00 



(15) 



Apart from a constant, Eq. ( 15 1 is a linear function of the 



creation and annihilation operators, hence, it is a linear 
function of the collective atomic quadrature variables 

X„ = ia,, + al)/V2 and P„ = (a„ - aJ,)/iV2. (16) 



Namely, Eq. (|T5| can be written as 
6 



7VOoo + \/2]V^ (ReO„o^„ + ImO„oP„). (17) 

a^O 

Let us choose real numbers Sind ^i, and the single- 
particle basis vector \(j>i) such that 



(18) 



Normalization implies that £md are the mean and 
variance of the single-particle operator, respectively. 



e, = {c^o\d^'^'\c^o)-eo^iAo^X- 



Eq. (17 1 now simphfies to O 



(19) 
(20) 



We 



can assign to O a single oscillator quadrature variable 



0-{0)o 



(21) 



and in the reference product state of the system, by con- 
struction, this collective atomic oscillator is in the ground 
state and (AXi)g — \. 

As an example, we mention that the natural normal- 
ization factor for the quadrature operators Q, that are 
assigned to the transverse angular momentum compo- 
nents, is not the spin quantum number J nor the macro- 
scopic expectation value of Jx- Rather, it is related to the 
variance of Jy and in the reference state, respectively, 



J I, 



p. 



Jz 



y2(AX)I' 



(22) 



Here (AJy)o = (AJ2)q — J /2 only for the coherent spin 
state, but not for a generic reference state. As long as 
the reference state is a product state and there is no en- 
tanglement among the particles, the corresponding col- 
lective atomic oscillators are in the ground state. Finally, 
we note that the Heisenberg uncertainty relation implies 
[Xy^P^J\/i < 1. The reason why Xy and may have a 
non-canonical commutation relation is discussed in the 
next subsection. 



C. Operators acting on different oscillators 

Let us now consider two collective operators, A and B, 
that can be linearized in the above manner. As we have 
shown, we can assign to them the oscillator quadrature 
operators Xa and Pb, respectively, such that 



A^{A)o+j2iAA)lX 



B=(B) 



2{AB)IPb. 



(23) 
(24) 



Note that the position or momentum quadratures of the 
oscillator modes are fully equivalent and the distinction 
depends on the choice of the single-particle basis and, in 
particular, the phase of the basis vectors. For the two 
operators, it remains a question, however, whether we 
can choose the basis in such a way that Xa and Pb are 
quadratures conjugate to each other acting on the same 
oscillator, or whether they belong to two independent 
atomic oscillators. To answer this question in general, 
we introduce the unnormalized vectors 



(B«-(i?W)o)|0o), 



(25) 
(26) 
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and analyze the following three cases. 

a. Parallel case. If \a) and |6) are parallel to each 
other (l^*) = Ae"''|a)), then Xa and Pb belong to the same 
atomic oscillator. Indeed, we can set = |a)/||a||, so 
that 



Pb = Xi cos ip + Pi sin ip, 



(27) 



where tp — arg(a|6). 

For example, in the vicinity of the fully polarized 
spin coherent state (|0o) = \F,mp = F)), the trans- 
verse components of the total angular momentum oper- 
ator, J = X^fLi F^^\ define the quadratures of the same 



atomic oscillator. The two vectors in Eqs (251 and (26 1 



|2/)sFW|0o) = 072|i^,m| 
|^)^fW|0o)=*|2/) 



= F-1) 



(28) 
(29) 



are parallel to each other with cp — 7r/2, so Xy — Xi and 
= A, and 



(30) 



In Appendix Ja 



we show in general that Xa and Pb 
defined in Eqs ( 23 1 and ( 24 1 are conjugate quadratures of 
the same atomic oscillator if and only if |0o) is a minimum 
uncertainty state with respect to ^'^^ and B^^\ 

b. Orthogonal case. If \a) and \b) are orthogonal 
to each other {{a\b) — 0), then Xa and Pb belong to 
completely different atomic oscillators. Choosing |0i) = 
|a)/||a|l and |</)2) = -i\b)/\\b\\ yields 



X, 



Xl: 



Pn 



P2 



(31) 



In recent experiments [D[lji the uncertainty in was 
significantly reduced by coherently squeezing the spin of 
each individual spin-4 cesium atom via two-axis counter- 
twisting [3] described by the Hamiltonian, 



The expectation value (T(^^)o is zero, since it is zero in 
the initial x-polarized spin state, and T'^^ is a constant 
of motion in the process of internal squeezing. It is easy 
to see that {Fz^^)o is also zero. Noting that T^-^^ changes 
the azimuthal quantum number of the hyperfine sublevel 
by two, whereas Pi^^ changes the same quantum number 
by only one, we immediately see that |t) is orthogonal 
to \z). Therefore, completely different atomic oscillators 
are accessed by the transverse components of the total 
angular momentum and the twisting operator. 

When the reference state is the a;-polarized coherent 
spin state and we take the standard spin projection eigen- 
basis \(j)a) — \F, mf, — F — a), the oscillators associated 
with Jz and T are the first and second ones, respectively. 
The idea of using both these oscillators in a quantum 
memory and a proposal on how to experimentally access 
the latter by means of stimulated four-photon processes 
is presented in Rcf. 4 . 

We remark here that any process that acts coherently 
on each atom (e.g., coupling to a classical field) changes 
the reference state itself but does not influence the atomic 
oscillators around the reference state. In this sense, inter- 
nal spin squeezing as a coherent process does not squeeze 
any of the collective atomic oscillators. In particular, 
it does not squeeze the quadrature operator Pz that is 
assigned to the z component of the angular momen- 
tum. It does reduce the noise in J^, but only because 
Jz « [2N{AfP)1]^^^Pz and the normalization factor in 
the brackets is reduced in the new reference state. In 
Sec. |III| we will investigate spin squeezing in more de- 
tails. 

c. General case. If \a) and |6) arc linearly indepen- 
dent but not orthogonal to each other, then A and B 
act on different but not independent oscillators. We can 
define an orthogonal basis via Gram-Schmidt orthogo- 
nalization and we can thus define Xa = Xi and 



Pb = {Xi cos p + Pi sin (p^ cos d + P2 sin i?, 



(36) 



0)2 



F. 



N 

E 



{P(^),P«}. (32) 



Thus each atom is prepared in the same state 



fe=0 



Ck\F,m^p^F-2k), 



(33) 



where ip — arg(a|&) as previously and cos?? — 
|(a|6)|/(||a||||6||) describes how parallel the two vectors 
are. The quadrature operators in this case have a non- 
canonical commutation relation, [Xa, Pb]/''' < 1- See 
Appendix [A| for a derivation of the formulae. 



D. Motion in rotating frame 



with only even number of spin excitations. In such a 
reference state, the quadrature operators assigned to the 
twisting operator T and to the transverse angular mo- 
mentum Jz belong to different atomic oscillators. Fol- 
lowing the prescription above, we introduce the corre- 
sponding vectors 



ii)^(f(i)_(TW)o)i0;;'^)-TW|0^<^ 



(34) 
(35) 



In a typical experiment, the atomic spins are placed 
in a homogeneous magnetic field where they precess co- 
herently and independently of each other. The reference 
state then also precesses, and we will describe the collec- 
tive excitations as deviations from this time dependent 
state. 

In general, we can split the Hamiltonian into two parts, 
H ~ Hq + Hi, one that acts coherently on each parti- 
cle and another that describes, e.g., interaction with an 
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external quantum field. Using the interaction free single- 
particle evolution operator, UQ^\t) = Texp—i Jq Hq^\ 
we can introduce the time dependent single-particle ba- 
sis \(l)a{t)) — UQ^\t)\(j)a) that defines a set of rotating 
atomic oscillators. Then we can define time dependent 
occupation number states analogously to Eq. (|3| , and the 
creation and annihilation operators ([s]) are all rotating 
accordingly in the Schrodinger picture, for example. 



(37) 



The rotating occupation number states satisfy the inter- 
action free Schrodinger equation. In the absence of Hi, 
the state of the rotating atomic oscillators and the num- 
ber of excitations remain unchanged. 

In the Heisenbcrg picture, the free Hamiltonian dis- 
appears from the equation of motion of the rotating 
atomic oscillator operators. For example, aanit) = 
U^t)aait)U{t) satisfies 



dt 



(38) 



To solve Eq. ( 38 1 , we need to express Hi in terms of the 



rotating oscillator variables. We assume that the inter- 
action does not bring the system far from the rotating 



reference state |^o(*))- Then Eqs (14 1 and (15 1 can 



be used to express the symmetric collective atomic opera- 
tors in Hi. Consider, for example, that Hi is an arbitrary 
function of the collective operator O. Then we write 

O « NOooit) + ^ [O„o(i)al(i) + H.c] 



E 



[0„/3(i) - Ooo{t)Sc.0]aUt)ap{t) (39) 



in place of O, where the matrix elements are 

OaMi) = {cl^M'^Ht)0^'^U^'\t)m. (40) 

Depending on whether \(j)oit)) is an eigenstate of O'^-* 
or not, we may neglect the quadratic terms in Eq. (39 1. 



Often Hi is at most quadratic in the oscillator operators, 
and the equation of motion reduces to a set of ordinary 
linear differential equations that can be solved. 



III. SPIN SQUEEZING 

For ensembles of spin-i particles, spin squeezing nec- 
essarily involves correlation among the spins ^5j. When 
these particles are actually the fictitious spins constitut- 
ing the real spin-.F atoms (e.g., the valence electron and 
the nucleons), correlations may be both intra- and in- 
teratomic. In the former case, the atoms are internally 
squeezed independently of each other, while in the latter 
case, there is genuine multi-atomic correlation. In this 
section, we investigate how the two ways of squeezing 
are related to each other. 



A. Internal spin squeezing 

We start our analysis with the situation in which the 
quantum uncertainty in the transverse component of each 
atomic spin, Fz''\ is reduced independently. No entangle- 
ment is created among the atoms in this way. To exem- 
plify our analysis, we will consider an ensemble of cesium 
atoms in the hyperfine level F = 4 of the atomic ground 
state (iSi/2 that is illuminated by an intense laser field. 
The light-atom interaction is off-resonantly tuned to the 
6>S'i/2-6P3/2 transition, hence photons are not absorbed, 
and transitions among the hyperfine sublevels as well 
as the ground-state energy shifts are mediated only by 
transfer of photons between different polarization com- 
ponents. Namely, for the jth atom interacting with a 
light field propagating in the x direction, the tensor light 
shift reads 

^^1=. 900(0) + .9il.(0)Fi^") 
+ .92 {-m^'^' + •^-(O)Fif + s+{0)F^'}') , (41) 

where the effective coupling constants gi sum up all the 
possible contributions from the different hyperfine levels 
of the relevant excited states and depend on the detuning 
as well. The quantum mechanical Stokes vector compo- 
nents describing the polarization state of light are 

•S+ = + iSy = ajjcii, 



(a^ai? - a[aL)/2, 



(42) 



where cir and are annihilation operators for right 
and left circular polarized, spatially localized pho- 
tons, and they obey the standard commutation relation 
[ai{z),aj{z')] = cSijS{z - z'). 

The photon flux commutes with any other operator 
in Eq. (41) and is usually treated as a c-number. The 



quadratic Fx^^"^ component can be compensated with the 
second order Zeeman shift. 



Hf=uj^F^^+fiF^^^ 



(43) 



by tuning the magnetic field such that = 52'/> [Illl]- To 
meet the two-photon resonance condition, the right and 
left polarized light components have to oscillate on the 
two Larmor sidebands. After summing over all atoms and 
performing the rotating wave transformation, we have 

H = .gis,(0)4 + 252 [sMy + h{^)T\ , (44) 
with the Hermitian collective operators ([32l and 



N ^ N 



We start from the a;-polarized spin coherent state and 
take the single-particle angular momentum eigenbasis 
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FIG. 1: Internal spin squeezing, equivalent to squeezing of a 
single F — 4 spin. The squeezing parameters xo, Co, and 
are shown as function of the integrated interaction strength, 
K — g2 J cj>{t) dt. The squeezing parameters are, in general, 
defined as follows: ~ '2{^JzY / {N F) is the uncertainty of 
the transverse spin component with respect to that in the orig- 
inal CSS [3], C'^ = 2(A Jz)^/{ Ja;) is the uncertainty compared 
to that of a CSS with the same longitudinal mean spin [5], 
5^ = 2NF{AJzf/{J^f is the noise in spin angle [7]. 



F 



a 



The operator T in Eq. (44 1 is 



responsible for squeezing F^-'^ (and anti-squeezing Fy-'^) 
of each atom, and we can select it by setting Sy — cj) /2 to 
have macroscopic expectation value, while the other two 
light operators, and s^, can be completely neglected. 
This corresponds to equally strong right and left circu- 
larly polarized light components, i.e., a strong linearly 
polarized field (rotating at twice the Larmor frequency 
in the laboratory frame). 

Under the action of Hq = g2<t>T, the initial spin coher- 
ent state evolves into the product state \Ui^\t)\4>ci)] 
in which every atom is internally squeezed, while the 
atomic oscillators expanded around this rotating refer- 
ence state are in their vacuum states. Therefore, the sta- 
tistical properties of collective observables are completely 
determined by the single-particle expectation values. In 
particular, the noise in the transverse component of the 
total angular momentum is AJ^ = \/~N/!^Fi^\ where 
Ai^i^-* is the uncertainty in the transverse spin compo- 
nent of a single atom, and squeezing has the same limit as 
for a single spin-i^ particle. The final degree of squeezing 
depends only on the dimensionless time integrated inter- 
action strength K = g2 J 4>{t) dt, that is proportional to 
the total number of photons in the squeezer pulse. Fig. [l] 
shows the degree of squeezing as function of this param- 
eter. 



B. Projection based spin squeezing 

In the next step, we consider a collective process in 
which the noise in Jz is reduced via quantum non- 
demolition (QND) measurement and feedback [5J [HI UHl 
HU m m [11 [T5J Hg. This kind of squeezing creates 
interatomic entanglement. To measure Jz in a non- 



destructive way, we shall couple it to a "meter" sys- 
tem, e.g., a light field propagating along the z-axis. This 
probe field is far detuned, so that we can neglect the sec- 
ond order light shift. We also neglect spontaneous emis- 
sion, absorption of the probe beam, and other sources 
of imperfections that may actually limit spin squeezing 
[m [T5] m [TS] . The interaction Hamiltonian in this con- 
figuration is Hi = giSz{0)Jz- 

The Stokes vector component Sz is a QND variable 
in this interaction: it is not modified by the interaction 
itself while the light passes different segments of the sam- 
ple. It is therefore conventional to treat the accumulated 
interaction as if the ensemble of atoms as a whole were 
interacting with a single light mode integrated along the 
pulse, Si{t) = Jq"^ Si{ct — ^,t)d^. For a strong classi- 
cal amplitude populating the a:-polarized light compo- 
nent, the meter system is the y-polarized quantum field 
and its quadrature operators can be defined as Xl = 
Sy/{cy/Np/2) and Pl = Sz/{c^J Np/2), where Np is the 
total photon number in the probe pulse. Initially, the me- 
ter is in the vacuum state and (AXij')^ = (AP[")2 = i. 

Regarding the atomic ensemble, we start from the 
generic reference product state |0o) ^^^d assume that 
the system stays in the vicinity of this state. We further 
assume that the polarization in this reference state points 
in the x direction {{Jy)o — {Jz)o = 0) and that |(/)o) is not 
an eigenstate of J; 



II B 



According to our results in Sec 
we assign to Jz a collecti ve atomic oscillator with quadra- 
ture variable A = Jz/ \/2{AJz)l. This is the atomic os- 
cillator the light is coupled to and, for the moment, it 
is enough to consider only this mode. It is initially in 
the vacuum state with (AXj")^ = (APi")^ = |. In this 

generalized Holstein-Primakoff approximation. Pi is also 
a QND vari able. The effective interaction Hamiltonian 
is Hi = gi^/Np{AJz)lPiPL/T, where T is the duration 
of the probe pulse, and the following input-output rela- 
tion holds for the quadrature operators in the Heisenberg 
picture. 



^out 

^out J^^^ 



HP] 



L J 



pout ^11] 



kP" 



r>out 



= P, 



L 1 



(46) 
(47) 



where k = gi ^ Np{AJz)f^ is the time integrated coupling 
that quantifies the strength of the measurement. The 
meter system is then read out by measuring Con- 
ditioned on the measurement outcome x, the atomic os- 
cillator becomes squeezed. The new means and variances 
are [15] 



(AXr*)2 = (1 + k2)/2^ (48) 



(AP°"*)2 = 



1 1 



21 



The uncertainty in Jz thus becomes 



(49) 



(50) 
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To ensure that the ensemble is still polarized in the x 
direction, the measurement result is fed back. Ideally, the 
feedback consists of a momentum displacement generated 
by the position quadrature Xi . In the x-polarized coher- 
ent spin state or whenever |(/)o) is a minimum uncertainty 
state with respect to Fy^^ and Fz^^ (such as the internally 



squeezed states in Sec. Ill A), the transverse components 



of the total angular momentum are orthogonal quadra- 
tures of the same atomic oscillator, irrespective of the 
degeneracy of the atomic ground state. Namely, Jy is 
proportional to Xi, so any Hamiltonian proportional to 
Jy (e.g., a magnetic field applied in the y-axis or a circu- 
larly polarized light beam propagating in the y direction) 
suffices to accommodate the ideal feedback. 

For a generic \4>q), however, a different feedback opera- 
tion is required. The quadrature operator assigned to Jy 
is now a combination of Xi and Pi, and there is also a 
contribution from an independent second atomic oscilla- 
tor, as described by Eq. @. Both (J™') and (J°"*) are 
proportional to the measurement outcome, that is, the 
transverse component of the total angular momentum 
is slightly tilted in the yz plane. The feedback appara- 
tus should be appropriately reoriented to correctly cancel 
this transverse angular momentum component. We also 
note that the mean spin ( J^,) may get further reduced if 
the feedback procedure excites a second atomic oscillator. 



C. Combining the two ways of spin squeezing 



b) 





t 



FIG. 2: (color online) a) Imaginary setup for combining in- 
ternal squeezing with measurement and feedback. The pulse 
propagating in the x direction is responsible for internal 
squeezing, the other one is the probe that measures Jz. b) 
Pulse diagram of the two kinds of squeezing applied in se- 
quence (1,2) and simultaneously (3). 



We have seen that spin squeezing can be achieved ei- 
ther by squeezing the corresponding atomic oscillator or 
by coherently acting on the reference product state of 
the system. Here we discuss how the two methods can 
be combined. Let us imagine the following oversimpli- 
fied scenario (Fig. |2|. The atomic ensemble, that is ini- 
tially prepared in the x-polarized coherent spin state, is 
illuminated (sequentially or simultaneously) by two light 
pulses: (i) The squeezer pulse, propagating in the x di- 
rection (i.e., in the direction of the atomic polarization), 
realizes the interaction Hamiltonian Hq = g24'T and 
gives rise to internal squeezing as detailed in Sec. |III A 



The relevant parameter of this pulse is the time inte- 
grated interaction strength K = g2 J <j>{t) dt. (ii) The 
probe pulse, propagating in the z direction, has a pho- 
ton flux (j>p{t). The QND interaction with the atomic 
sample is given by Hi = giSz(0)Jz- The Stokes vec- 
tor component Sy{z) of the outgoing field is continuously 
measured as in Sec. |IIIB| The relevant parameter of 
the probe pulse will be the effective integrated coupling 
= ^gfNF J (j)p{t) dt. Given the two pulses, we analyze 
the following three combinations: (1) internal squeezing 
is followed by measurement based squeezing, (2) the same 
sequential squeezing but in reverse order, and (3) when 
the two methods are applied simultaneously. 



1. Internal squeezing followed by measurement 

Let us first address the case in which an internally 
squeezed ensemble is further squeezed by QND measure- 
ment. In Eq. (50 1, that already applies for a generic refer- 
ence state, (AJz)q = -^FN for the a;-polarized coherent 
spin state, but {AJz)q < \FN for internally squeezed 
states. Therefore, better squeezing can be achieved if 
the measurement is preceded by internal squeezing. How- 
ever, the strength of the measurement also depends on 
the reference state. In the same measurement setup (i.e., 
same probe pulse length and intensity), the time inte- 
grated coupling is reduced if the ensemble is internally 
squeezed: k = hxoi where k = <7i NpFN /2 denotes 
the coupling for the x-polarized coherent spin state, and 
Xo = \/'^{^Jz)f)/{NF) is the squeezing parameter of the 
internally squeezed reference state. The effective inte- 
grated coupling k does not depend on the actual state 
of the atomic ensemble, it characterizes only the probe 
pulse. The final degree of squeezing, parametrized by 



2{AJzy/{NF), then reads 



xl{K,k)^[Xo'{K) + K' 



(51) 
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FIG. 3: Combination of the two ways of spin squeezing, 
as function of the integrated interaction strength of internal 
squeezing. Curve xi shows the final degree of squeezing when 
internal squeezing is followed by measurement and feedback, 
while curve X2 corresponds to the reverse order. Applying 
the two processes simultaneously results in curve xs- The ef- 
fective time integrated coupling of the measurement is k = 2 
for all the three cases. As a reference, the dashed line (xo) 
shows internal squeezing only. In the inset: evolution of the 
non- vanishing matrix elements Im J^o for a = 1, 3, 5, and 7 
(with decreasing line thickness). 



as function of the integrated interaction strength of the 
initial internal squeezing process, K, and the effective in- 
tegrated coupling of the measurement, k. As a compari- 
son, the curve xo in Fig. [3] shows internal squeezing only. 
The difference between xi and xo is due to the second 
phase of squeezing, namely, to the QND measurement. 
The contribution of the measurement to the overall de- 
gree of squeezing decreases with internal squeezing, since 
the coupling k = kxo also decreases. 



2. Measurement and feedback followed by internal squeezing 



Let us now consider the two squeezing processes in the 
reverse order. We start from the x-polarized coherent 
spin state, and consider the single-particle angular mo- 
mentum eigenbasis |0q,) = nip = F — a) as the initial 
basis for the atomic oscillators. In this basis, Eq. (30) 
holds and Pi is the quadrature assigned to Jz- This 
quadrature is first squeezed by measurement and feed- 
back, so (APi)^ becomes as given by Eq. (49 1, with 
K = k. Then we switch off the probe pulse and switch on 



the squeezer pulse. Since the Hamiltonian Hq — g24>T 
realizes a coherent action on each atom, the state of the 
collective atomic oscillators are not changed, only the un- 
derlying basis is rotated, \4>a{t)) = Uo{t)\(j)a)- In this new 
basis, however, no longer corresponds to the squeezed 
quadrature Pi. Instead, we have 



/2]V^ImJ^oW^'a- 



(52) 



The coefficients J^o(t) = {(j)a{t)\Fi^\(t)o{t)) depend on 
time only through the integrated interaction strength K. 



In what follows, we will write K instead of t in the ar- 
gument. The functions J^q{K) are shown in the inset of 
Fig. [3j Since the atomic oscillators are independent of 
each other in this basis, the variance of Jz is given by 



{AJzY = 2NJ2 [Im J^oW]'(AP^)' 



(53) 



Taking (APi")2 = 1(1 + ^2)-! according to Eq. (|49|, and 
(AP™)^ = i for a > 1, the final squeezing parameter is 



xl{K,k)^xliK) 



[ImJfo(if)]V2 



(54) 



We can see in Fig.[3]that the effect of the first, projection 
based squeezing decreases with the duration of the sec- 
ond, internal squeezing, and the minimum of the overall 
squeezing is reached before the minimum of xo- 



3. Measuring while squeezing internally 

Finally we analyze the case when the two kinds of in- 
teraction are applied simultaneously. We assume that 
the two pulses are switched on at the same time and, for 
simplicity, we will consider a constant intensity profile 
for both pulses. The integrated interaction strength of 
internal squeezing, K, as well as the effective integrated 
coupling of the measurement, k, can then be controlled 
via the intensities of the respective pulses. 

Since the atomic oscillator basis rotates due to internal 
squeezing, the probe pulse couples at different time in- 
stants to different combinations of the collective atomic 
oscillators, as specified in Eq. (52 1. Therefore, it is not 



enough to consider a single atomic oscillator and a single 
light mode integrated along the pulse. In Appendix |B] 
we formulate the proper Maxwell-Bloch equations of mo- 
tion and derive a differential equation for the covariance 
matrix of the collective atomic quadratures. Given the 
time evolution of the reference frame, we can numerically 
evaluate the squeezing parameter (see xs in Fig. [3]) . 

Finally, we note that the presence of an external mag- 
netic field, that is essential for realizing the internal 
squeezing Hamiltonian, actually prohibits the QND mea- 
surement of Jz ■ This is due to the fact that the transverse 
angular momentum precesses in the yz plane. This prob- 
lem can be circumvented by using two oppositely oriented 
cells of atoms as in 17j. 

To conclude this section, we have shown that the two 
kinds of squeezing reduce each others efficiency, and their 
effect do not simply add up. If internal squeezing is ap- 
plied first, the coupling strength of the QND readout is 
decreased with respect to that in a coherent spin state. 
If, on the other hand, the measurement based squeezing 
is applied first, internal squeezing will mix the collec- 
tive atomic oscillators, and it is no longer the originally 
squeezed oscillator that corresponds to the angular mo- 
mentum component. 
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IV. APPLICATION TO QND MEMORIES 

In atomic quantum memories based on ofF-resonant 
QND interaction [17l [HI [HI EQI [21], the ensemble is 
usually prepared in a coherent spin state, and the atoms 
are effectively treated as spin-i particles, for which the 
conventional Holstein-Primakoff approximation leads to 
a straightforward oscillator description. We now gener- 
aHze this approach to arbitrary atomic level structures, 
and we investigate the interplay between the choice of 
reference internal state and the collective variables. 

Let us now briefly recall the functioning of the quan- 
tum memory based on QND measurement and feed- 
back [17 , but with emphasis on more general reference 



states. We consider the same QND setup as in Sec. IIIB 



the atomic ensemble, that has a mean spin pointing to 
the X direction but is otherwise prepared in a generic ref- 
erence state, is illuminated by a strong x-polarized light 
field propagating in the z direction. The photonic quan- 
tum oscillator, whose state we wish to map on the en- 
semble, is the copropagating y-polarized quantum field 
integrated along the pulse. Irrespective of what the ref- 
erence state actually is, the information is stored as a dis- 
turbance around the arbitrary reference state, and this 
disturbance can be described as excitation of the collec- 
tive atomic oscillator which is assigned to Jz in the given 
reference state. 

After the light has first passed the sample, the input- 
output relation among the quadrature variables of this 
collective atomic oscillator and those of the photonic os- 
cillator is readily given by Eqs (|46]) and (|47l. The X^"* 
quadrature of the outgoing light field is then measured, 
and the measurement outcome is fed back. The atomic 
variables after the feedback become [TFl [55] 



X'l = Xi" -I- KPf , Pi" 



1 



X 



L ' 



(55) 



and both X™ and PJf have now been transferred to the 
collective atomic oscillator variables Pi and Xi, respec- 
tively. The transfer of information is, however, not per- 
fect: -q^ = 2(AX{")^/k^ units of vacuum noise is added 
to the original P™ quadrature due to the initial uncer- 
tainty in the collective atomic quadrature X™. It has 
been proposed to mitigate this imperfection by squeez- 
ing the atomic spin state before the memory operation. 
Such squeezing operation could be performed, for exam- 
ple, by an additional QND measurement of Jy. Indeed, 
this would create interatomic correlations and squeeze 
the collective atomic oscillator. It has also been sug- 
gested |2] that internal squeezing can enhance the fidelity 
of the quantum memory as well. We have shown, how- 
ever, that internal squeezing does not reduce AXi, so it 
is not cleax whether it really enhances the performance 
of the memory. 

We have seen in Sec. [Ill C 1| that the integrated cou- 
pling K = KXzo decreases when is squeezed internally, 
while it increases if Jz is anti-squeezed, for example, due 



to internal squeezing of Jy instead of Jz- Assuming for 
a moment that |0o) is a minimum uncertainty state and 
(AJy)o(AJ2)o = NF/A, we obtain for the amount of ad- 
ditional noise 



7/ 



2(AXi")2 



Xy 



(56) 



where Xy is the squeezing parameter for Jy and k is the 
effective coupling of the QND mapping. We thus observe 
that the fidelity of the quantum memory may improve if 
Jy is internally squeezed. This does not happen because 
the collective atomic oscillator is squeezed, but rather be- 
cause the reference state |</)o) is a minimum uncertainty 
state, and the adjoint, anti-squeezed, Jz implies an in- 
creased coupling K. 

We remark here that (A Jz)§ and, thus, the coupling k 
can be increased even without squeezing Jy. The refer- 
ence state is, however, not a minimum uncertainty state 
in such a case, so the transverse angular momentum 
components no longer correspond to the same collective 
atomic oscillator. This fact has two consequences. The 
first regards the feedback and has already been addressed 
at the end of Sec. [IIIB] by reorienting the feedback appa- 
ratus, the measurement outcome can be eliminated from 
the state of the first collective atomic oscillator, but at 
the cost of exciting a second one. The second conse- 
quence concerns the read-out of the memory. 

The memory can be read out by a similar procedure, 
interchanging the role of the atomic and photonic oscil- 
lators. When a read-out light pulse passes the atomic 
ensemble, the quadrature operators are transformed into 



X'i 
X 



R 



X'l 



Pi - Pi. 



P" 



pu 



(57) 
(58) 



Then a measurement of the atomic Xi quadrature and 
a subsequent feedback onto the light pulse should fol- 
low. In Ref. jlH], an additional light beam was suggested 
to couple a "meter" system to the angular momentum 
component Jy . If the reference state is not a minimum 
uncertainty state, however, we may not have access to 
the Xi quadrature through Jy. In a generic reference 
state. 



jy cx (Xi cos (p + Pi sin (p) cos ^ + X2 sin d, 



(59) 



with tan 95 = -(pi^^)o/(T(i))o and cos^ = ((Pi^')^ + 
(T(i))§)/[4(APj$^^)2(APi^')2]. The A component can 
be ruled out by reorienting the measurement device in 
the same way as for the feedback. The contribution from 
the second atomic oscillator X2, however, cannot be elim- 
inated. Instead of measuring Xi^ the best we can actu- 
ally measure is x' = X[" + jH^^f. The outgoing Ught 

quadrature then turns into P^^' = jj{x' - X'{ - j^X^') > 
and the quadratures of the read-out pulse after the feed- 



10 



back on the light read 



"V^OUt 



X 



R ' 



rjout 
^R 



K 



P'r' 



tan -d 
KCOS ip 



-X" 



(60) 
(61) 



Besides the noise term in Eq. ( 56 ) and the noise intro- 



duced by the Xfi quadrature of the readout pulse, an 
additional noise term appears if |(/)o) is not a minimum 
uncertainty state. Assuming that both the atomic os- 
cillators are initially in their vacuum states, the amount 
of noise in the P'j^* quadrature (in vacuum noise units) 
reads 



V 



(1)\2 



> 



1 



(62) 



where the lower bound is obtained using the Cauchy- 
Schwarz inequality ( A5 1 . 



To summarize this section, we have analyzed the oper- 
ation of the quantum memory based on QND interaction 
and feedback in a generic reference state. We have shown 
that, although internal squeezing does not squeeze the 
collective atomic oscillators, it may enhance the perfor- 
mance of the memory by enhancing the coupling strength 
of the QND interaction. We have also pointed out the 
difficulties arising when the reference state is not a min- 
imum uncertainty state with respect to the transverse 
angular momentum components. 



V. CONCLUSIONS 

If each spin-i particle in an ensemble is prepared in 
the same single-particle state, it is always a spin coher- 
ent state, and small perturbations to the product state 
can be well described by a single collective oscillator de- 
gree of freedom in the Holstein-PrimakofF approximation. 
In this paper, we have introduced a generalization of this 
method to describe an ensemble of d-level atoms in the 
vicinity of an arbitrary product state (not necessarily 
spin coherent state) . We have defined d — 1 independent 
collective atomic oscillator modes, and we have specified 
how to express collective operators (namely, permuta- 
tion invariant sums of single-particle operators) in terms 
of the oscillator creation and annihilation operators. 

We have applied our formalism in particular to spin 
squeezing of atoms. We have analyzed two different 
methods: internal squeezing and QND measurement 
based squeezing, and identified the collective oscillators 
that have become squeezed. We have shown that the 
two kinds of squeezing reduce the effect of each other. 
When the atoms are first internally squeezed, the cou- 
pling strength of the QND measurement is reduced and, 
therefore, the measurement based squeezing is not so ef- 
ficient. If we first project the ensemble's state and then 
continue with internal squeezing, then the latter process 
will mix the collective atomic oscillators and, at the end 



of the process, it is no longer the originally squeezed oscil- 
lator that corresponds to the transverse angular momen- 
tum component. We have also considered the case when 
the two kinds of squeezing is applied simultaneously. 

Finally, we have analyzed a quantum memory scheme 
for storing quantum states of light in atomic ensembles 
based on QND interaction and feedback using an arbi- 
trary initial product state. We have shown that internal 
squeezing can reduce the noise of the memory as an in- 
direct effect because of an enhanced coupling strength, 
caused in fact by the anti-squeezing of the adjoint atomic 
spin component. 



APPENDIX A 



In this appendix, we derive the formulae in Sec. II C 
We define the single-particle basis |0i) and 102) from the 
vectors \a) and \b) using Gram-Schmidt orthogonaliza- 
tion, |0i) = |a)/||a|| and 



. \b)-{{a\b)/\\ar)\a) 
V\\b\\'-\{a\b)\y\\a\\ 



(Al) 



The matrix elements in Eq. (|17| are 

^10 = l|a|| 

^20 = 0, 



Sio = (a|6)/||a||, 

B,o = ^V\\br-M 



Vll«ll 



(A2) 



Introducing the mixing angles ip = arg(a|6) and cosd 
|(a|6)|/(||a||||6||), we can write 



Bw = ||fe||e''^cosi?, B2o = i\\b\\smi}. 



(A3) 



Given that ^2(AB)§ = V2N\\b\\, after substituting 
Eq. (A3 1 into ([T7| and comparing it to Eq. (24), we arrive 
at Eq. (36 1, which we wanted to prove. 



Now we show that it is sufficient and necessary for the 
quadrature operators Xa and Pb to belong to the same 
atomic oscillator and to be conjugate to each other that 
|(/)o) is a minimum uncertainty state with respect to the 
single-particle operators A^^") and B'^^K From Eq. (36 1 



and from the definition of the mixing angles, we have the 
commutation relation 



[Xa,Pb]/i = 



lm{a\b) _ i([iW,B(i)])o 



|a||||6|| (AA(i))o(AB(i))o" 



(A4) 



From the Heisenberg uncertainty relation we know that 
the absolute value of the real number at the right-hand 
side of Eq. ( A4 1 is less than or equal to 1, and the inequal- 



ity is saturated, by definition, for minimum uncertainty 
states. Exactly for such states = ±i, and the 

two quadratures are canonically conjugate to each other. 
The Cauchy-Schwarz inequality. 



|a||||6|| > x/(Re(a|6))2 + (Im(a|6))2, 



(A5) 



then implies that \a) and \b) are parallel to each other 
{cosd = 0) and {a\b) is pure imaginary (ip — if and 
only if \(f)Q) is a minimum uncertainty state. 
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APPENDIX B 

Here we derive the degree of squeezing when measure- 
ment based and internal squeezing is simultaneously ap- 
plied. See Fig. |2] and Sec. |III C| for a description the 
setup. 

Internal squeezing, governed by the Hamiltonian Hq — 
g24>T , acts coherently on each atom. The corre- 
sponding time dependent single-particle basis \4>a{t)) = 
ij^\t)\(l3a) defines the rotating atomic oscillators. The 
coupling to the meter system is described by the in- 
teraction Hamiltonian Hi = giSz{'0)Jzj where Sz now 
refers to the probe pulse, whose photon flux is (j)p{z, t) = 
(j)p{t — z/c). With the coordinate change £_ = ct — z, it 
is convenient to consider the propagating slices of the 
y-polarized quantum field of the meter system. The 
quadrature operators of the slice, that enters the sam- 
ple at the time instance ^/c, are 



^ Sy{ct - ^,t) 

^ v/<^p(e/c)/2' 



^ Szjct ~ ^,t) 
~ ^m/c)/2-' 



(Bl) 



in the Heisenberg picture. The equation of motion for 
the Heisenberg operators of the rotating atomic oscilla- 
tors is given by Eq. ( 38 ) . First we express the inter- 



action Hamiltonian with the quadrature operators using 
Eqs ([52| and (IbTI), 



Hi it) « 51 iV0p it) J2 Ini Jao W (t) Pl {ct, t) . (B2) 
Then we arrive at the following Maxwell-Bloch equations: 



dt 



P„(t) =0, 



dt 



dt 



Xo.{t) - giJNcbpit)lmJ^^{t)pL{ct,t), 



L (e, i) = .91 ^JNMt) E "^"0 it) it) edict -0. 

To analyze the effect of the continuous light measure- 
ment on the atomic oscillators, let us divide the probe 
pulse into short segments of duration r. The quadrature 
operators of the segment, that enters the sample at the 
time instance ^/c, are X^{t) = J^^'^'^ XLiC jt) d^ /ic^/r), 

and similarly for P^{t). We assume that the change in 
the reference state due to internal squeezing, as well as 
the change in the probe photon flux (j)p{t) can be ne- 
glected during the passage of a single light segment, and 
that the evolution of the atomic system can be obtained 
by sequential interaction with the segments. Right after 
the segment has passed the sample, the following input- 
output relation holds. 



pout 
^out 



P"an £)Out 
a 1 -ff 



Pi' 



(B3) 
(B4) 
(B5) 



where the labels "in" and "out" respectively mean be- 
fore and after the passage of the light segment in 
consideration, and = gi^J 4>pi^/ c)tN is an effec- 
tive coupling constant. Introducing the vector y = 
A, • ■ • , X2F, P2F, P^)^, we can write Eqs (|B3])- 



(B5) as a matrix equation y — Sy™. For Gaussian 



states, as in the case considered here, the system is fully 
characterized by the vector of expectation values (y) and 
the covariance matrix Fy = ((y^ — {yi))iyj — iVj))}, that 
transform as (y°^') = S(y'") and T™* = Sr'"S^, respec- 
tively. 

When the light segment enters the sample, it is com- 
pletely uncorrelated with the atomic ensemble: the 
incoming covariance matrix is block diagonal, F'" = 
diag(r"J, 1), where Fat is the covariance matrix of the 
collective atomic oscillators and 1 is the 2x2 identity 
matrix describing the initial vacuum state of the light 
segment. We denote the outgoing covariance matrix by 



pout Q 

qT pout 



(B6) 



where C is a 2 x 2P matrix describing the light-atom 
correlations. After the interaction, the meter system is 
measured. Conditioned on the measurement outcome, 
the effect of the light measurement on the atomic covari- 
ance matrix is given by the relation j23j 



r' 



1 



pout] 



1 




(B7) 



Combining the formulas above, we can write a differ- 
ence equation that describes the change in the atomic co- 
variance matrix due to the weak QND measurement |15j . 
If the segment is short enough, it is sufficient to keep 
only the leading order in and in the limit of in- 
finitesimal r, we arrive at a set of ordinary differential 
equations for the atomic covariance matrix T. In our 
case, we have a closed subset of equations for the matrix 
elements describing momentum-momentum correlations, 
7„^^((P„-(P„))(P^-(P^))), 



dt 



7a/3 



^E7aa'ImJ^,oIm>/|'o7/3'/3, (B8) 

a'P' 



where k^/t = gfN(j)p. If the time dependent coefficients 



Im J^Q(t) are known, Eq. (B8) can be solved. Then we 
can calculate the uncertainty in the transverse angular 
momentum component ( 52 ) and the squeezing parameter 



(B9) 



a/3 



Fig. [3] shows the final degree of squeezing obtained by 
numerically integrating Eq. ( B8 1 for constant intensity 



profiles for both the squeezer pulse and the probe pulse 
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